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ABSTRACT 

This  Investigation  analyzed  families  of  ascent 
trajectories  to  determine  total  velocity  require- 
ments, division  of  velocity  impulses  and  the  ex- 
pected errors  in  position  and  velocity  upon  arrival 
in  the  terminal  region  caused  by  uncertainties  in 
the  position,  magnitude  and  direction  of  the 
velocity  vector  at  burnout  of  the  booster  stages. 

Conditions  at  burnout  of  the  booster  and  apogee 
of  the  trajectory  were  solved  using  an  IBM  650 
computer;  the  error  analysis  was  accomplished  by 
numerical  integration  on  IBM  709  equipment  repeated 
up  to  the  apogee  time  of  the  zero-error  trajectory. 

From  the  tabulated  and  graphical  results  pre- 
sented herein,  some  trade-offs  between  propulsion 
requirements  and  the  accuracy  of  guidance  and  con- 
trol systems  can  be  seen. 
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OBJECT 

It  is  the  object  of  this  thesis  to 
present  in  graphical  form  the  effects 
of  errors  in  position  and  velocity  at 
booster  burnout  on  the  position  and 
velocity  of  a  rendezvous  vehicle  at 
zero-error  apogee  time. 


CHAPTER  I 


INTRODUCTION 


Requirements 

Long  before  the  launching  of  the  first  earth 
satellite  the  military,  scientific,  economic  and 
political  ramifications  of  such  a  vehicle  were  con- 
sidered and  discussed.  Each  of  the  above  considera- 
tions generates  adequate  reasons  for  a  satellite 
rendezvous . 

The  most  important  and  probably  the  most 
immediate  requirement  is  in  the  field  of  the  mili- 
tary. In  orbit  about  the  earth  today  are  satellites 
capable  of  reconnaissance;  little  additional  equip- 
ment is  required  to  extend  this  capability  to  ferret 
and  jamming  missions  or  ICBM  early  warning.   Inter- 
rogation and  possible  negation  of  such  enemy  vehicles 
are  missions  that  require  a  satellite  rendezvous. 

Satellite  stations  in  space  require  the  physical 
exchange  of  personnel,  equipment  and  fuel.  Com- 
pletion of  such  an  exchange  requires  a  "soft"  rendezvous 
and  continuation  of  the  two  vehicles  in  the  same  orbit 
until  the  transfer  is  complete. 


3 
Both  civil  and  military  •rganizations  are 
developing  communications  satellites;  these  de- 
vices involve  much  complex  and  expensive  equip- 
ment that  will,  in  time,  require  repair  or  re- 
placement.  Satellite  rendezvous  is  a  solution 
to  this  problem . 

Communications  satellites  are  given  special 
consideration  in  this  report  by  the  inclusion  of 
data  for  orbit  altitudes  as  great  as  20,000  nauti- 
cal miles  above  the  surface  of  the  earth.  Earth 
satellites  much  above  this  altitude  are  of  little 
interest  since  they  are  subject  to  large  pertur- 
bations by  the  sun  and  moon  and  hold  no  advantage 
over  the  satellite  with  a  period  of  24  hours.  The 
effect  of  such  perturbations  and  the  difficulties 
to  be  encountered  in  maintaining  a  high  altitude 
orbit  are  discussed  by  Edelbaum  in  reference  1. 

Background 

The  subject  of  satellite  rendezvous  has  re- 
ceived its  share  of  attention  in  the  past  several 
years.  Literally  hundreds  of  papers  and  articles 
have  been  written  on  the  various  phases  of  rendez- 
vous and  interception.  Most  of  these  works  are 
very  general  and  deal  primarily  with  ideas  rather 
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than  figures  and  hardware.   Other  authors  look  into 
the  problem  in  detail  and  find  it  to  be  enormous; 
it  seems  to  take  a  well-funded  project  to  generate 
any  amount  of  quantitative  date  on  a  particular 
phase.  Many  studies  are  initiated  on  the  basis 
that  the  proposed  payload  and  vehicle  character- 
istics are  known  and  the  task  is  to  find  out  how 
(or  if)  the  desired  mission  can  be  accomplished. 
Much  of  this  work  is  security  classified  and  is, 
therefore,  unaccessible.  It  is  the  purpose  of 
this  thesis  to  conduct  an  academic  yet  practical 
study  of  satellite  rendezvous  so  that  students  may 
gain  an  insight  into  a  very  complex  problem. 

An  important  part  of  this  subject  that  is 
specifically  not  included  in  this  thesis  is  the 
topic  of  guidance  and  control  during  the  terminal 
phase  of  the  rendezvous  maneuver.  This  topic 
appears  to  be  a  study  in  itself  inasmuch  as  the 
terminal  phase  shifts  to  "cloifed  loop"  control 
with  feedback  from  radar  as  the  rendezvous  vehicle 
closes  on  the  target.  Some  facts  about  the  capa- 
bilities of  the  acquisition  system  and,  in  fact, 
entire  specification  of  the  means  to  be  used  by 
the  rendezvous  vehicle  in  completing  the  terminal 
phase  must  be  outlined  to  justify  the  definition 
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of  a  successful  arrival.  Sears  and  Pelleman  in  re- 
ference 2  give  an  excellent  example  of  this  type 
study.   Nason  establishes  the  conceptual  limitation 
of  proportional  navigation  in  reference  3. 


CHAPTER  II 


STATEMENT  OF  THE  PROBLEM 


General 

The  earth  satellite  is  a  unique  device 
since,  unlike  vehicles  that  operate  in  the  earth's 
atmosphere,  its  flight  is  a  regular  pattern  in 
space  determined  primarily  by  the  velocity  imparted 
to  the  vehicle  at  the  beginning  of  its  orbit.  While 
velocities  of  earth  satellites  are  high  compared  to 
any  motion  within  the  earth's  atmosphere,  the  motion 
of  the  satellite  is  predicted  with  accuracy  and  it 
is  this  statement  that  forms  an  important  base  uti- 
lized in  this  thesis. 

Satellite  rendezvous  is  the  problem  of  bringing 
together  two  objects  of  relatively  small  size  when 
the  original  distance  between  the  two  objects  and 
their  velocities  are  both  very  large  by  ordinary 
engineering  standards.  When  the  target  orbit  ele- 
ments are  known,  the  task  of  the  rendezvous  vehicle 
(RV)  is  to  match  them.  In  practically  any  circum- 
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stances,  even  the  simplest  case,  there  exist  an 
infinite  number  of  possibilities  for  accomplish- 
ing a  rendezvous . 

The  general  case  of  matching  the  ephemeris 
of  an  arbitrary  target  after  launch  from  an 
arbitrary  point  in  space  (surface  of  the  earth) 
is  unmanageable  and  probably  unrealistic .  Al- 
though an  effort  is  made  to  keep  the  problem  as 
general  as  possible,  certain  conditions  are  set 
before  proceeding . 

The  Specific  Problem 

The  study  begins  with  the  assumption  that  the 
orbit  of  the  target  and  the  trajectory  of  the  RV 
are  in  the  same  plane .  If  the  orbit  of  the  target 
has  an  inclination  less  than  the  latitude  of  the 
launching  point  then  plane -change  is  required  as 
the  target  never  passes  over  the  launch  point . 
With  angles  of  inclination  greater  than  this  lati- 
tude the  launch  point  is  in  the  plane  of  the  target 
tvfice  each  day  giving  two  opportunities  for  launch 
into  the  target  plane.  Except  for  polar  orbits, 
one  of  these  opportunities  requires  a  westward 
launch,  an  undesirable  choice,  thus  the  desirable 
eastward  launch  into  the  target  plane  is  available 
only  once  a 
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Now  Consider  further  the  odds  that  the  target 
is  at  a  suitable  position  for  rendezvous  at  the 
time  the  launch  point  passes  through  the  target 
orbit.  Obviously  a  wait  of  many  days  may  be  re- 
quired before  the  correct  conditions  are  met . 

A  study  of  the  wait  time  problem  is  made  by 
Houbolt  in  reference  4.  In  this  paper  he  de- 
monstrates for  a  specific  case  the  limited  number 
of  rendezvous  opportunities  available  in  a  given 
period  and  points  out  the  increase  in  this  number 
given  by  the  "Adjacent  Launching  Technique".  In 
this  scheme  the  RV  is  launched  into  an  orbit  of 
the  same  altitude  as  the  target  with  the  target 
at  a  suitable  position  for  rendezvous.  When  the 
RV  crosses  the  target  orbit  it  adapts  to  this 
new  orbit  with  a  plane-change  maneuver  and  begins 
the  terminal  phase  of  the  rendezvous. 

Inasmuch  as  changes  of  orbit  inclination  are 
accomplished  by  an  impulse  normal  to  the  original 
velocity  vector,  the  change  of  plane  requirement 
is  in  addition  to  the  in-plane  requirements  and 
is  not  considered  further. 

Since  the  target  is,  by  definition,  the  point 
in  space  to  which  we  wish  to  move  the  RV,  it  is 
reasonable  to  assume  that  the  target  has  a  regular 
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motion  easily  predicted.  A  circular  orbit  best 
meets  this  requirement  and  this  is  not  an  un- 
realistic assumption;  for  many  reasons  an  approxi- 
mately circular  orbit  is  desired  for  reconnaissance, 
survey  or  communications  satellites. 

Now  within  a  given  plane  and  with  a  circular 
target  orbit  there  still  remain  an  unlimited 
number  of  trajectories  that  meet  the  requirement 
of  placing  the  RV  at  the  satellite's  position  in 
space.  These  trajectories  form  the  basic  conic 
sections,  ellipse,  parabola  and  hyperbola,  passing 
through  the  two  required  points.   Only  a  minimum 
of  computation  is  required  to  see  that  hyperbolic 
and  parabolic  trajectories  have  excessive  velocity 
requirements  when  compared  with  the  elliptical 
path.  The  elliptical  trajectories  that  qualify  are 
those  with  apogee  greater  than  or  equal  to  the  al- 
titude of  the  target.  The  former,  often  called  a 
"fast  ellipse",  crosses  the  target  orbit  and  this 
is  undesirable  for  two  reasons s 

1.  To  match  the  orbit  of  the  target  the  RV 
requires  in-plane  normal  thrusting  at  the  cross- 
over point.  Normal  thrusting  is  inefficient  and 
is  only  used  as  a  last  resort. 
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2.  Closing  velocities  at  the  cross -over  are 
very  high  and  rendezvous  requires  critical  timing. 
The  ascent  ellipses  considered  in  this  thesis 
are  those  with  apogee  equal  to  the  altitude  of  the 
target.  This  type  trajectory  gives  a  flight  path 
at  apogee  that  is  tangent  to  the  orbit  of  the 
target . 

Within  the  conditions  of  being  in-plane  and 
tangent  at  arrival/  there  is  an  alternative,  four- 
impulse,  "parking  orbit"  ascent.  Briefly,  parking 
orbit  ascents  consist  of  a  launch  into  a  low  alti- 
tude orbit,  a  waiting  period  until  the  relative 
positions  of  RV  and  target  are  correct,  then  a 
transfer  to  the  target  orbit.  The  primary  advantage 
of  such  a  scheme  is  the  relaxation  of  strict  re- 
quirements for  leaving  the  launching  pad  at  a  given 
time;  the  low  altitude  orbit  allows  the  RV  to  "gain" 
on  the  target.  The  disadvantages  of  such  a  system 
are  related  to  this  thesis  in  that  errors  in  burnout 
velocity  are  not  only  multiplied  by  the  requirement 
for  two  "transfer"  maneuvers,  but  initial  errors  in 
many  cases  "grow"  during  the  wait  period.  These 
and  other  factors  are  considered  in  detail  by  Duke, 
Goldberg  and  Pfeffer  in  reference  5. 
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Coast  ing  in  free  flight  from  burnout  of  the 
booster  stages  until  the  KV  is  at  the  altitude 
of  the  target  is  termed  a  "direct  ascent".   Such 
a  trajectory  has  a  minimum  velocity  requirement 
in  the  classical  Hohmann  maneuver  and  an  upper 
limit  in  the  vertical  ascent.  This  thesis  will 
be  concerned  with  these  two-impulse  type  rendez- 
vous trajectories. 

To  review,  then,  the  specific  problem  in- 
vestigated in  this  thesis  is  an  in-plane  rendezvous 
with  a  target  in  circular  orbit  utilizing  a  direct 
ascent  elliptical  trajectory  that  is  tangent  to  the 
target  orbit. 
Basic  Assumptions 

The  earth  is  assumed  to  be  a  homogeneous 
sphere.  The  earth  is,  of  course,  quite  irregular; 
it  is  more  an  ellipsoid  than  a  sphere. 

This  ob lateness  has  three  main  perturbating 
effects  upon  the  orbit  of  an  earth  satellite,  all 
dependent  upon  altitude  and  inclination  of  the 
orbit : 

1.  Regression  of  the  nodes.  This  is  motion 
of  the  orbital  plane  about  the  earth's  axis  in  a 
direction  opposite  to  that  of  the  satellite. 
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2.  Precession  of  the  line  of  apsides.  This 
motion  is  an  angular  rotation  of  the  axis  of  the 
orbit  in  inertial  space.  It  may  be  in  the  same 
or  the  opposite  direction  of  the  satellite  motion 
depending  upon  the  amount  of  inclination. 

5.  The  oblateness  of  the  earth  also  causes 
slight  fluctuations  of  the  radius  (altitude)  and 
period  of  the  satellite. 

All  of  the  above  effects  are  important  in 
sustained  satellite  operations  and  their  effects 
on  orbits  adds  to  the  knowledge  of  the  shape  of 
the  earth,  however,  such  effects  are  of  small 
magnitude  in  direct  ascent  rendezvous  trajectories. 

As  an  example,  in  reference  6  the  formula: 

A  -  2» (1.657  x  10"-*)  = 75     COS  1 

(Re  ♦  z8)2 

is  given  as  the  measure  of  angular  motion  of  the 
nodes  in  radians  per  revolution  of  the  satellite. 
For  an  altitude  of  1000  n.  miles  and  with  a  cir- 
cular orbit  inclined  at  an  arbitrary  angle  of  45°, 
the  regression  of  the  nodes  is  4.38  x  10"*  radians 
for  every  revolution  of  the  target.  The  ascent 
trajectories  described  in  this  thesis  are  of  short 
duration  compared  with  the  period  of  the  satellite 
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target  and  only  a  portion  of  this  regression  of 
the  nodes  takes  place  while  the  RV  is  in  flight. 
For  the  1000  n.  mile  target  the  ascent  time  varies 
from  5134  seconds  to  843  seconds  or  11#  to  40#  of 
the  period  of  the  target.  Thus  the  regression  of 
the  nodes  is  from  13,000  to  47,000  feet  at  this 
geocentric  radius. 

This  example  is  given  only  to  show  the  order 
of  magnitude  of  the  motion  in  this  case.  In  the 
rendezvous  situation  the  perturbations  affect 
the  RV  as  well  as  the  target  at  the  same  time  and 
in  the  same  direction.  The  error  is,  then,  of 
differential  order  when  compared  with  the  above. 

If  these  effects  are  expressed  as  functions 
of  position  they  can  be  introduced  into  the  inte- 
gration of  the  equations  of  motion  for  the 
rendezvous  computations, but  they  are  not  the  sub- 
ject of  the  present  investigation  and  can  influence 
it  only  through  higher  order  effects. 

To  continue  with  the  basic  assumptions, 
booster  burnout  altitude  is  taken  as  sixty  nautical 
miles  for  all  cases.  Without  a  detailed  examination 
of  the  booster  performance  and  maneuvering  scheme 
it  cannot  be  determined  whether  this  altitude  is 
appropriate  for  a  given  trajectory, but  some  arbi- 
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trary  altitude  is  needed  for  a  common  base  from 
which  to  compare  the  coasting  performance.  This 
altitude  is  assumed  in  numerous  references  such 
as  Eggleston  in  reference  7. 

It  is  further  assumed  that  the  ephemeris  of 
the  target  is  known  exactly.  The  technique  of 
tracking  an  earth  satellite  and  determining-  its 
ephemeris  is  fairly  well  formalized  and  high 
speed  digital  computer  routines  are  available  for 
this  purpose.  The  technique  is  described  by  Sandorff 
in  reference  8.  For  this  study,  however,  uncertainties 
of  the  satellite  position  are  simply  uncertainties  in 
the  origin  of  the  final  error  coordinates .  While  this 
may  be  an  Important  consideration,  it  is  independent 
and  can  be  introduced  as  necessary  in  a  system  study. 

Although  the  ability  to  predict  launch  time 
continually  improves  with  experience  in  missile 
operations,  there  are  still  uncertainties  in  time 
of  lift-off  as  well  as  error  build-up  during  burn- 
ing. The  time  of  flight  or  coasting  time  varies 
considerably  with  the  particular  trajectory 
selected  (as  seen  below)  giving  some  flexibility 
in  matching  time.  A  pre-launch  computer  keeps 
track  of  this  time  and,  over  some  useful  range 
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determined  by  other  considerations,  selects  the 
exact  trajectory  to  match  time.   In  a  similar 
manner,  during  the  burning  period  performance  is 
monitored  and  cut-off  corrections  computed  to 
match  the  time.  In  very  general  terms,  Duke, 
Goldberg  and  Pfeffer  in  reference  5  take  notice 
of  this  type  error  but  its  implementation  de- 
pends upon  system  specifications.  If  an  error 
in  time  alone  is  accepted  and  no  compensation 
is  made  then  the  arrival  error  is  easily  com- 
puted from  the  satellite  velocity.  To  get  to 
the  point  of  the  study,  it  is  assumed  further 
that  booster  burnout  is  compelled  to  occur  ex- 
actly on  time. 
Uncertainties  Considered 

The  uncertainties  considered  in  this  thesis 
are  those  in  the  booster  burnout  velocity  vector. 
These  errors  are  described  as: 

1.  Position  errors 

a.  ahead  (RANGE  error) 

b.  above  (ALTITUDE  error) 

c.  to  one  side  (TRACK  error) 

2.  Velocity  errors 

a.   error  in  magnitude 
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b.  error  in  direction 

(1)  deflected  above  or  below 
the  ideal 

(2)  deflected  to  one  side  of 
the  ideal 

The  exact  meaning  of  "position  errors"  de- 
pends upon  the  type  of  guidance  equipment  utilized 
by  the  booster  stage;  the  full  implication  of  this 
point  is  carefully  discussed  in  Chapter  III. 
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CHAPTER  III 

COMPUTATIONAL  METHODS  AND  PROCEDURES 

Analytical  Computation  of  Reference  Trajectories 
The  elements  that  describe  the  target  orbit 
and  the  family  of  possible  direct  ascent  ellipses 
tangent  to  it  are  defined  with  the  help  of  Fig.  1 
below. 


Non-iit\©A     WltXtff 


Circular   McVtUVfo 
Velocv^ 


Fig.  1.  Definition  of  Elements  Used  in  Computational 
Procedures . 
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Prom  the  figure  it  is  seen  that  selecting 
a  target  altitude  fixes  R^  for  all  ascent  ellipses, 
since  the  rendezvous  condition  is  that  apogee 
occurs  at  target  altitude.  An  additional  variable 
is  needed  to  specify  a  particular  trajectory  from 
the  possible  family  of  ellipses, and  the  range 
angle,  4>  ,  best  describes  what  is  a  long  or  short 
ascent.  Although  these  two  variables  are  selected 
for  classification  of  results,  the  conventional 
two-body  equations  are  more  familiar  in  terms  of 
the  standard  parameters  e,  a  and  9.  The  desired 
values  are  computed  by  successive  substitution  in- 
to the  numbered  equations  below. 

These  equations  are  not  involved,  but  they 
require  careful  numerical  evaluation  for  accurate 
results.  With  dimensions  in  millions  of  feet, 
the  IBM  650  carrying  eight  significant  digits  is 
adequate  for  computation  directly  in  engineering 
units.  The  three  physical  constants  for  the 
numerical  work  are: 

G   (gravitational  constant  times  mass 
m  of  the  earth) 

1.40768  x  10l6ft5/sec2    Ref.  8 

R   (radius  of  a  sphere  with  volume 
equivalent  to  that  of  the  earth) 

5457o866  nautical  miles     Ref.  9 
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International  nautical  mile  -  £076.1155  feet 

Ref.  10 
Altitude  traversed  in  ascent 

z  ,   »  z   -  z  , 

<*      si 

Minimum  eccentricity 


(3.1) 


*2L 


2^  +  z* 


min 
Angle  from  perigee 


9  -  (180°  -  <*>  ) 


(3.2) 


(3.3) 


Eccentricity 

cosO  +1            ,\ 

1_ 
4 

2 
cos9,  +  1 

COS9+1+  =— 

L                      Rl 

C     ■                         "" "                    "  ""'!_"    '          T  I 

2  [    COS9+1+  ^J      \ 

Satellite  velocity 


V. 


Gm 


Rl  +  z* 


Booster  burnout  velocity 


V 


s 


+  1-e 


'*+A 


cos9+l+  g- 

s.  1  J 


(3.4) 


(3.5) 


(3.6) 


Apogee  velocity 


V. 


Vsf^ 


(3.7) 
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Injection  velocity  increment 

V2  -  Vs  -  V*  (5.8) 

Semi -major  axis 


*1  +  ZoC 
a  "  1  +e 

Booster  burnout  attitude  angle 


(5.9) 


Time  to  apogee 

a  -  R-,  r  a  -  R 


a 


Yfe  {T  "  arccos(     ae  1)  +  e  sin[arccos(     ae  1) 

(5.H) 

This  information,  computed  for  12  target  alti- 
tudes from  200  to  20,000  nautical  miles  and  range 
angles  of  180,  90,  45  and  22.5  degrees  permits 
selection  of  a  smaller  representative  sample  for 
detailed  error  analysis .  Target  altitudes  of  400, 
1,000,  5,000  and  20,000  n.  miles  are  used  in  the 
error  analysis  to  follow. 

The  Error  Analysis 

A  complete  error  analysis  of  even  this  basic 
set  of  trajectories  by  analytical  methods,  if 
possible  at  all,  is  certainly  complex  because 
several  orbital  elements  are  often  changed  by  a 
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single  error  at  booster  burnout.  For  instance, 
the  variation  in  semi -major  axis  can  be  found  as 
a  function  of  a  change  in  )f  ,  but  before  an  error 
can  be  isolated,  the  corresponding  change  in  ec- 
centricity and  the  shift  of  the  line  of  apsides 
must  be  determined.  Finally,  arbitrary  position 
as  a  function  of  time  is  trancendental. 

Braham  and  Skidmore  in  Ref .  11  give  analytical 
expressions  for  changes  in  altitude,  final  velocity 
and  final  flight  path  angle  as  a  function  of  certain 
initial  condition  errors.  These  expressions  are 
cumbersome .  The  one  for  final  flight  path  angle 
requires  evaluating  six  lengthy  expressions  to  get 
the  necessary  partial  derivatives,  which  in  turn 
require  evaluating  four  other  equations  equally 
as  long  to  solve  for  some  of  the  coefficients. 
Even  if  a  complete  set  were  available,  evaluation 
in  any  quantity  calls  for  a  machine  computer  pro- 
gram. 

Such  analytical  solutions  are  interesting 
from  the  theoretical  point  of  view,  but  it  is 
felt  that  the  most  practical  way  to  get  a  com- 
plete error  analysis  is  by  numerical  integration 
of  the  equations  of  motion. 
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The  motion  of  a  body  subject  to  an  accelera- 
tion is  described  by  the  vector  equations: 


V  -   J  a  dt  +  VQ  (3.12) 

R  -   J  V  dt  +  RQ  (3.13) 


In  coasting  flight  near  the  earth,  the  only 
acceleration  is  that  due  to  the  gravitational 
field,  G,  which  is  a  function  of  R.  Thus,  for 
step-by-step  integration,  the  two  simultaneous 
equations  are  actually  solved  in  order  as  follows: 


Gm 
,3 


(3.14) 


R' 

This  acceleration  is  expressed  in  components  in 
whatever  coordinates  V  is  given,  and  each  com- 
ponent of  V  and  R  is  computed  as: 

Vi-G1  At  +Vl0  (3.15) 

Ri  "  vi  At  +  Rio  (3.16) 

The  resulting  R  is  used  to  compute  a  new  G;  the 
resulting  V  is  the  V0  for  the  next  step,  and  the 
resulting  R  is  also  the  R  for  the  next  step.  This 
procedure  is  repeated  up  to  apogee  time. 

A  complete  program  with  this  capability,  de- 
veloped by  Mr.  Robert  Brown  of  the  M.I.T.  Instru- 
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mentation  Laboratory  and  adapted  to  this  problem, 
is  used  in  the  error  analysis.  With  this  program 
and  the  IBM  709  computer  at  the  M.I.T.  Computation 
Center,  trajectories  for  four  altitudes  and  four 
range  angles  are  thoroughly  investigated. 

The  method  is  straightforward  and  accurate 
when  used  with  a  suitable  integrating  time  step. 
The  validity  can  be  checked  for  each  trajectory 
by  integrating  over  the  complete  time  of  flight 
with  no  errors  introduced  and  comparing  the 
apogee  conditions  with  the  analytical  results 
from  the  equations  above.  As  an  example  of  the 
accuracy  attained  the  zero -error  numerical  in- 
tegration for  ascent  to  a  5,000  n.  mile  altitude 
at  a  range  angle  of  90°  using  a  50  second  time 
step  agreed  to  within  80  feet  ALTITUDE  and  7  feet 
RANGE  distance  -  this  is  at  a  geocentric  radius 
of  over  51  million  feeto 

The  program  uses  spherical  coordinates  in 
the  actual  integration,  but  keeps  track  of  the 
components  in  terms  of  an  orthogonal  frame  called 
a  platform  frame.   This  frame  is  introduced  to 
simplify  the  presentation  of  the  error  components, 
and  has  a  positive  X  axis  in  the  RANGE  direction 
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along  the  horizontal  In  the  direction  of  flight, 
a  positive  Z  axis  upward  along  the  vertical 
ALTITUDE  direction,  and  a  positive  Y  axis  to  the 
right  looking  in  the  direction  of  flight. 

A  platform  frame  is  specified  because  for 
inertially  guided  vehicles  directions  are  re- 
ferred to  a  stabilized  platform.  The  platform 
is  non-rotating  in  earth-centered  Inertial  co- 
ordinates (type  l)  or  its  Z  axis  tracks  the  local 
vertical  (type  2).  The  effect  of  position  error 
on  both  types  is  shown  in  Fig.  2. 

The  first  type  seems  most  likely  for  this 
application,  but  it  is  used  here  for  a  more  basic 
reason.  With  error  effects  available  for  the 
first  type,  they  can  be  computed  for  the  second 
type.   Note  from  the  figure  that  ALTITUDE  errors 
alone  are  the  same  in  either  system,  but  that 
with  RANGE  or  TRACK  error,  the  relative  alignment 
of  the  G  vector  depends  on  the  type  system  assumed 
It  is  seen  though,  that  a  type  2  RANGE  error  con- 
sists of  a  type  1  RANGE  error  plus  a  velocity  mis- 
alignment (  SY  )  which  is  easily  computed  from 
booster  burnout  radius  and  the  RANGE  error. 

Effects  due  to  the  actual  change  in  height 
caused  by  a  type  1  RANGE  or  TRACK  error  are  negli- 
gible because,  for  reasonable  errors,  this  change 
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in  height  is  very  small,  and  the  corresponding 
position  error  derivatives  are  also  very  small. 
The  same  is  not  true  for  misalignment  errors. 
The  errors  may  be  fractions  of  a  milliradian, 
but  this  kind  of  error  has  a  large  end  effect. 
Finally,  for  TRACK  error,  the  motion  due 
to  the  type  1  position  error  must  be  obtained 
by  numerical  integration,  but  for  the  type  2 
platform  where  the  ALTITUDE  component  of  velocity 
is  aligned  with  the  local  vertical,  the  entire 
motion  lies  in  the  plane  of  the  error  trajectory. 
When  a  trajectory  lies  entirely  in  a  plane  and 
the  inclination  with  respect  to  another  plane 
is  known,  the  lateral  displacement  is  just  a 
function  of  the  geometry.  This  point  is  brought 
out  in  more  detail  in  Chapter  IV. 
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local 


Comparison  of  two  types  of  platform  coordinate 
frames.   Angle  for  reasonable  position  error 
greatly  exaggerated. 


JVW, 


Ref.  local  horizontal 


lv,'\  ■  IV.l 


Ref.  local  horizontal 


Error  in  magnitude  of 
booster  burnout  velocity 


In-plane  misalignment 
of  booster  burnout 
velocity. 


re<ev«nct 

local 
vertical 


Track  Error  -  Type  1 
system  as  used  herein. 


Lateral  Misalignment 
of  booster  burnout 
velocities . 


Fig.  2.   Platform  Coordinates  and  Definitions  of  Booster 
Burnout  Errors . 


27 


CHAPTER  IV 


RESULTS 


Summation  of  Velocities 

The  summation  of  booster  burnout  velocity  and 
injection  velocity  increment  (XV)  for  various  alti- 
tudes and  range  angles  is  shown  in  Fig.  2.  It  must 
be  borne  in  mind  that  the  values  represented  here 
are  the  theoretical  minima;  they  do  not  reflect  the 
additional  velocity  that  must  be  provided  to  over- 
come drag  and  gravity  losses  during  accelerations 
to  V±o 

What  is  considered  important  about  these  plots 
is  the  effect  of  altitude  and  range  angle  on  XV. 
Note  that  while  analysis  always  prescribes  the 
Hohmann  maneuver  for  minimum  XV,  this  is  not 
necessarily  the  primary  consideration  in  a  practical 
case.   For  instance,  with  a  400  n.  mile  orbit  alti- 
tude the  range  angle  is  reduced  to  55°  before  XV 
increases  \%   over  the  Hohmann  value;  it  is  reduced 
to  18°  before  XV  increases  5#.  Similarly  for  a 
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1000  n.  mile  orbit  altitude  the  range  angle  is 
down  to  40°  before  XV  increases  by  5#. 

These  examples  are  cited  only  to  call  at- 
tention to  what  may  be  one  corner  of  a  many 
sided  compromise. 

Velocity  Division 

Pig.  4  is  a  plot  of  the  separate  components 
of  XV  versus  range  angle  for  a  number  of  orbit 
altitudes.  The  information  is  displayed  on  two 
separate  plots  for  low  and  high  altitudes  to 
eliminate  the  clutter  of  several  overlaps. 

V,,  the  booster  burnout  velocity,  is  the 
velocity  which  must  be  imparted  to  the  RV  to 
send  it  coasting  up  to  apogee.  V~,  the  injection 
velocity  increment,  represents  that  propulsive 
capability  that  must  be  carried  to  the  rendezvous 
point  in  order  that  satellite  velocity  be  matched. 

Prom  another  point  of  view,  V~  is  the  relative 
velocity  between  target  and  RV.   If  the  ascent 
trajectory  is  perfect  and  time  matched  exactly, 
the  injection  thrusting  can  be  initiated  by  elapsed 
time,  but  with  error  buildup  during  the  ascent  it 
seems  most  likely  that  the  terminal  phase  will  be 
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closed-loop.  Thus,  this  relative  velocity  may 
be  important  during  the  acquisition  oif  the  target 
by  the  RV. 

As  before,  for  computation  of  propulsive  re- 
quirements the  drag  and  gravitation  losses  must 
be  added  to  V-,°  Then,  once  an  effective  exhaust 
velocity  or  specific  impulse  is  known  or  assumed, 
quick  preliminary  calculations  can  be  made  con- 
cerning vehicle  design.  These  might  consist  of 
selecting  an  altitude  and/or  range  angle  within 
the  capability  of  a  certain  system,  or,  in  deter- 
mining what  configuration  and  staging  or  what  re- 
start capability  is  needed  to  accomplish  a  given 
mission. 

This  information  forms  the  second  part  of 
the  propulsive  picture,  and  another  side  of  the 
overall  compromise  involved. 

Error  Derivatives 

The  results  of  the  error  analysis  are  re- 
duced to  partial  derivatives  since  it  is  shown 
that  the  variation  of  a  final  component  with  re- 
spect to  a  change  in  the  initial  value  of  any 
other  component  is  linear  over  what  is  considered 
to  be  a  generous  range  of  variation.  This  direct 
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proportionality  is  verified  for  the  1000  mile 
orbit  for  speed  errors  up  to  60  feet  per 
second,  for  angular  deviations  of  the  velocity 
up  to  two  milliradians  and  for  position  errors 
up  to  12,000  feet  in  any  direction.  As  shown 
for  the  errors  in  booster  burnout  velocity, 
plotted  on  Pig.  5,  there  is  no  tendency  to  de- 
part from  linearity  at  60  feet  per  second.  Be- 
cause the  functions  being  considered  are  smooth 
and  continuous,  they  do  not  bend  abruptly  so  it 
is  considered  that  booster  burnout  velocity 
errors  up  to  100  feet  per  second  and  angular 
deviations  involving  components  of  that  size 
could  be  entertained  with  confidence.  It  does 
not  seem  likely  that  such  large  errors  are 
realistic,  although  what  is  tolerable  certainly 
depends  on  what  is  said  about  the  terminal 
system  capabilities.  In  any  case,  the  additional 
computer  time  required  to  investigate  by  this 
method  and  establish  the  absolute  limit  of  appli- 
cability is  not  justifiable  with  the  many  other 
items  being  studied  with  the  computer. 

The  types  of  booster  burnout  uncertainties 
described  above  are  classifiable  as  in-plane  and 
lateral  errors.  With  three  components  each  to 
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define  the  position  and  velocity  vectors,  six 
differentials  are  sufficient  to  describe  any 
possible  error  in  the  booster  burnout  condition. 
To  express  the  reference  apogee  or  final  value 
of  each  of  the  six  vector  components  in  terms 
of  each  of  the  six  initial  components  leads  to 
36  partial  derivatives.  There  is,  however, 
clearly  no  way  for  the  in-plane  effects  to  in- 
duce any  lateral  perturbations  and,  as  expected, 
the  lateral  errors  introduce  only  infinitesimal 
effects  on  the  final  in-plane  performance.  With 
this  independence,  the  number  of  error  derivatives 
is  reduced  to  20. 

While  it  seems  more  natural  to  describe  the 
initial  velocity  errors  in  terms  of  variations 
in  magnitude  and  direction,  the  final  velocity 
errors  are  given  in  terms  of  orthogonal  components. 
This  is  done  primarily  for  ease  of  combination  of 
components  as  demonstrated  later.  Also,  since  both 
the  target  and  the  reference  apogee  point  have  zero 
Y  and  Z  components  of  velocity,  the  computed  error 
components  are  the  actual  relative  velocity  in  the 
ALTITUDE  and  TRACK  directions  as  viewed  from  the 
satellite. 

The  X  component  of  the  error  velocity,  while 
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generally  quite  small  compared  to  either 
the  actual  or  relative  velocities  in  the  RANGE 
direction,  is  still  important.  It  is  a  part  of 
the  terminal  problem,  and  the  rendezvous  is  not 
complete  until  all  components  of  relative 
velocity  have  been  exactly  matched. 

The  derivatives  relating  the  initial  and 
final  variations  are  made  up  of  combinations 
of  the  following: 

1.  In-plane  -  Initial  Position  (x  and  z) 
and  velocity  (v,  and  "Y  ) 

Final  position  (X  and  Z)  and  velocity 
(Vx  and  Vz ) . 

2.  Lateral  -  Initial  Position  (y)  and  velocity  (f). 
Final  position  (Y)  and  velocity  (VY). 

The  general  form  of  the  drivative  as  used  for 
evaluating  ascent  errors  is ;$ initial ) '   This 
does  not  preclude  using  the  inverse.  If,  for 
instance,  it  is  desired  to  define  allowable  final 
error  first,  then  th©  permissible  booster  burnout 
deviations  of  a  type  could  be  found „  However, 
there  are,  in  general,  many  combinations  of  booster 
burnout  errors  that  can  combine  to  cause  similar 
final  errors. 
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The  derivatives  are  plotted  in  Figs.  6  through 
15  with  range  angle  as  the  primary  variable 
and  the  target  altitude  as  the  step  parameter. 
Each  is  explicitly  defined  on  the  respective 
plot  with  units  in  feet,  feet  per  second  and 
milliradianSo 

It  is  important  to  note  that  as  range 
angle  approaches  zero,  the  eccentricity  of  the 
trajectory  ellipse  approaches  unity  and  the 
regular  forms  of  the  equations  become  very 
sensitive.   Judging  from  a  few  of  the  cases 
for  which  the  vertical  ascent  may  be  evaluated, 
it  appears  that  the  error  functions  in  fact  are 
quite  unstable  at  low  range  angles.   For  this 
reason,  none  of  the  results  should  be  extra- 
polated below  the  22.5°  range  angle  where  the 
error  computations  were  cut  off.  There  seem 
to  be  practical  reasons  for  avoiding  very  steep 
ascents  such  as  the  high  relative  velocities, 
highest  possible  2.V,  and  need  for  rapid  change 
of  attitude  in  the  terminal  phase.   For  such 
reasons,  it  is  considered  that  the  lack  of 
coverage  at  these  short  angles  is  not  important. 
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Independence  of  Derivatives 

A  final  important  result  is  the  independence  of 
the  in-plane  derivatives.  To  investigate  all  possible 
combinations  of  booster  burnout  errors  in  varying 
magnitude  involves  endless  computation,  but  the  total 
error  is  shoWn  in  Table  I  to  be  made  up  of  the  sum  of 
the  errors  due  to  the  separate  effects.  Table  I  is 
compiled  by  introducing  a  random  combination  of  errors 
into  the  1000  n.  mile,  '90°  range  angle  trajectory  and 
checking  the  computed  output  for  component  addition 
from  the  separate  derivatives. 

TABLE  I 


Pinal  Error  Components 


Booster  Burnout 

X 

Z 

vx 

vz 

Error 

feet 

feet 

ft/sec 

ft/sec 

(1)  +20  ft/sec 

-44,295 

58,693 

-79.94 

40.25 

(2)  -2  mr 

-103,095 

40,404 

-73.15 

-10.03 

(3)  1  mile  ahead 

-7,711 

6,048 

-10.92 

-.75 

Sum  of  (1)  +  (2) 

-147,388 

99,097 

-153.09 

30.22 

Computed  Results 

-147,387 

98,929 

-153.10 

29.89 

Sum  of  (1)  +  (3) 

-52,004 

64,741 

-90.86 

39.50 

Computed  Results 

-51,994 

64,753 

-90.87 

39.48 

Sum  of 
(1)  +  (2)  +  0) 

-155,099 

105,145 

-164.01 

29.47 

Computed  Results 

-155,080 

104,964 

-164.04 

29.09 

48 

Let  the  assumed  booster  burnout  errors  be  de- 
signated by  gx,  Sz,  &V,  and  &Y  for  in-plane 
errors  and  by  by   and  Si*  for  lateral  errors.  Now 
if  the  final  errors  are  designated  by  AX,  AZ,  AV 
and  AVZ  for  in-plane  errors  and  by  AY  and  AVy  for 
lateral  errors  then  a  complete  final  error  for  a 
booster  burnout  error  in,  say,  speed,  range  and 
lateral  alignment  is: 


AX  - 

ax 

Sv, 

♦-ti** 

AZ   * 

Sv1 

+  -SS  *>* 

AV 

3Vx 
8V1 

Btx 

bv 

+^i^ 

AV8- 

9vz 

8v1 

AY  - 
AV  - 

s«r01 
4/ ST 

where  the  partial  derivatives  are  read  from  the 
graph  for  the  appropriate  altitude  and  range  angle. 
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Table  I  demonstrates  that  this  linearization 
introduces  no  more  than  1#  inaccuracy  in  the  final 
computed  values  for  a  practical  case.  Special 
notation  seems  unnecessary  since  the  procedure 
is  so  simple.   Any  combination  is  readily  written 
if  it  is  remembered  that  each  initial  in-plane 
error  causes  four  components  of  final  error,  two 
in  position  and  two  in  velocity.  The  lateral 
errors,  which  are  independent  of  the  in-plane 
error,  each  cause  two  final  errors,  one  in  TRACK 
position  and  one  in  TRACK  velocity. 

Lateral  Position  Error 

The  lateral  position  error  due  to  a  TRACK 
error  at  booster  burnout  rates  special  mention 
because  it  illustrates  an  effect  due  entirely 
to  the  definition  of  position  error  used,  that 
is,  the  type  of  platform  assumed.  The  easiest 
place  to  view  this  effect  is  at  a  range  angle 
of  90  degrees.  It  seems  at  first  that  for  a 
lateral  position  error  at  booster  burnout,  the 
line  of  nodes  between  the  plane  of  the  actual 
ascent  ellipse  and  the  reference  plane  would 
lie  90  degrees  away,  or  at  apogee  for  this 
case.  The  integrated  results  shown  in  Fig.  14 
show  this  is  not  the  case. 
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Pig.  16  below  illustrates  how  this  comes 
about  with  the  central  angle  due  to  a  reasonable 
position  error  greatly  exaggerated. 


Mote 


Reference,  plane. 


turnout    POsWion  .and 
Velocity  cornponervVs  art 
iri    plcx-V^orm  coordinates 


VRM1K  PoV\Vior\  error 


V,  M«*  ^/R~) 


Altitude,  component  ©f  V, 


■Y« 


X^>  \v\Vo  j».yer 


Fig.  160  Effect  of  Platform  Definition  on 
Lateral  Position  Error. 
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Note  that  although  no  lateral  velocity 
appears  in  the  platform  coordinates,  with  a 
platform  that  stays  aligned  with  the  reference 
vertical,  a  sidewise  component  of  velocity, 
V1  sin  "ft  (y/R_),  does  show  up  with  respect  to 
the  reference  orbital  plane. 
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CHAPTER  V 


CONCLUSIONS 


If  this  thesis  has  met  its  objective  then 
the  plots  showing  the  various  error  derivatives 
should  clearly  demonstrate  the  primary  con- 
clusions to  be  drawn  concerning  the  effect  of 
target  altitude  and  range  angle  on  these  de- 
rivatives . 

To  review  the  results  in  general: 

1.  For  the  low  altitude  rendezvous 
the  summation  of  booster  burnout 
velocity  and  injection  velocity  in- 
crement (ZV)  is  only  slightly  great- 
er than  the  Hohmann  value  until  the 
range  angle  is  well  below  90°.  The 
components  of  ^V  follow  this  same 
pattern. 

2.  The  in -plane  error  derivatives 
for  a  1000  n.  mile  target  altitude 
are  accurate  for  differentials  of 

+  60  ft/sec,  +  2  milliradians  and 
12,000  feet  in  any  direction. 
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3.  The  lateral  error  derivatives  are 
independent  of  the  in-plane  error  de- 
rivatives. 

4.  The  in-plane  partial  derivatives 
are  independent  of  each  other  and  there- 
fore the  total  error  components  can  be 
formed  by  the  rule  for  total  differentials 

5.  The  propulsive  requirements  for 
the  20,000  n.  mile  orbit  altitude  are 
shown  to  be  so  great  that  mult i -staging 
will  certainly  be  necessary. 

At  this  altitude  the  sensitivity 
to  error  is  so  great  that  special  guidance 
techniques  will  also  be  required. 
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